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A stable oscillon was found in the (classical) analysis of the SU(2) gauged Higgs 
model under a 2:1 ratio between the Higgs and W masses, both analytically and 
numerically. We extend these results by constructing oscillons for all values of the 
Higgs mass and examining their stability, finding them to be stable when the Higgs 
mass is almost equal to twice the W mass. We also investigate the dependence of 
the oscillon lifetime on the mass of the Higgs, seeing different behavior for different 
regions of the Higgs mass. Close to the 2:1 mass ratio the interaction of the two 
fields leads to either a stable oscillon or to its quick decay, based on its amplitude 
and the exact masses of the fields. Moving further away from this ratio the fields 
seize to interact and the oscillon disperses. Furthermore, we show that the oscillon 
lifetime scales quadratically with the small amplitude parameter which is linked to 
the small amplitude expansion that we use to construct the oscillon. Finally, we 
list all families of oscillon solutions that can be formally constructed in this model 
and discuss their properties. The results can be easily extended to other composite 
oscillons in particle physics and beyond. 

PACS numbers: Preprint MIT-CTP 4412 

I. INTRODUCTION 

The formation of long-lived localized packets of energy is a common feature in nature. 
The literature of mathematics and physics is full with cases where nonlinearities in field 
theories lead to stable localized configurations. Relativistic field theories have been known to 
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exhibit such phenomena We can divide them in time- independent and time-dependent. Their 
longevity can be due to conserved charges or an effect of the interplay between nonlinearity 
and dispersion. 

We consider a specific type of localized object called an oscillon. Oscillons are localized in 
space, but oscillate in time (hence the name) and they do not posses any known conserved 
charges. In general they are not completely stable, in fact analytical and numerical argu- 
ments show an exponentially small radiative tail [1-4]. However they are found numerically 
to have lifetimes that can be orders of magnitude longer than the natural timescales in the 
theory, such as the oscillation period. This can persist even when they are placed in an 
expanding background [5]. If one considers quantum mechanical effects, the radiation of the 
oscillons is increased [6], leading to a shorter lifetime. 

Most documented cases of oscillons involve toy model scalar field theories, where various 
analytical and numerical techniques have been used to describe their properties, such as 
their size and lifetime [7-9]. Perhaps the most physically interesting characteristic of oscil- 
lons is their "natural" emergence in a wide range of scenarios, including inflation [10-12], 
bubble collisions [13], an expanding system with quasi-thermal initial conditions [14] and 
the interaction of a double well system with a heat bath [15]. Finally composite oscillons 
have been found to arise in a hybrid inflation setting [16]. 

Recent attempts to find oscillons in the Standard Model lead to the investigation of the 
electroweak sector. An oscillon in the SU(2) gauged Higgs model was found both numerically 
[17] and analytically [18] when the Higgs and W fields have a 2 : 1 mass ratio. Also, numerical 
studies revealed a stable oscillon in the full SU(2) x U(l) sector [19] and [20]. Furthermore, 
simulations have shown that during the electroweak symmetry breaking in an expanding 
universe, oscillons can emerge and capture a significant amount of the energy of the system 
[21]. 

In [18] Farhi et al construct an oscillon in a restricted version of the SU (2) gauged Higgs 
model and calculate the oscillon lifetime as a function of the Higgs mass. We explain the 
lifetime behavior observed by Farhi et al. and thoroughly investigate the whole region of 
Higgs masses, constructing all possible oscillon solutions, in the small amplitude expansion, 
and characterizing their stability. Wherever possible we provide intuitive arguments, sup- 
ported by analytical calculations and numerical simulations. Our analysis can be applied 
to any case of composite oscillons, meaning that the use of the SU (2) model is in no way 
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restrictive! In fact reducing the analysis of a complicated system, such as the SU(2) model, 
to simple arguments makes our treatment easily transferable. 

Our analysis of of the oscillon's behavior for all values of the Higgs mass reveals regions 
of qualitatively different behavior. We formulate the solutions using a small amplitude 
expansion with parameter e < 1. The size of the oscillon is 0(1/ e) and the dynamics 
happens on a timescale of 0(1 /e 2 ). We formally construct new oscillons solutions, who are 
shown to be unstable. Close to the 2 : 1 mass ratio, a formal oscillon solution exists and is 
found to be stable for a region of Higgs masses, where the Higgs mass is parametrized as 
m 2 H = Am 2 v +e 2 Am 2 H . Outside the stable region, the two fields that form the oscillon interact 
and depending on the available decay phase space, the Higgs field transfers its energy into 
the W or vice versa. As we move further away from the 2 : 1 mass ratio, the oscillon is 
described as a sum of two non-interacting fields, both evolving dominated by dispersion. 

The rest of the paper is organized as follows. In Section II and III we present the model 
and formulate the small amplitude expansion in the various regions of the Higgs mass. The 
region of masses is divided using the small parameter e, characterizing each region by the 
deviation of the mass ratio between the Higgs and W masses from 2. In Section IV the 
oscillon profiles for each oscillon solution are numerically derived. In Sections V and VI, we 
take a closer look at the evolution of the oscillon, including its lifetime and the interaction 
between the two constituent fields, for all regions of the Higgs mass and the expansion 
parameter. Conclusions follow in Section VIII. 

II. MODEL 

We follow the analysis presented in [17, 18] and will repeat the basic results for complete- 
ness. We consider the Lagrangian of an SU(2) gauge theory coupled to a Higgs doublet in 
3+1 dimensions 
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The theory contains three vector bosons (W's) with mass mw — gv/2 and a Higgs particle 
with mass m# = v \/2\. 

In order to be able to search for non-trivial solutions to the equations of motion we will 
constrain the Lagrangian by using the spherical symmetry [22, 23]. The spherical Ansatz is 
given by expressing the gauge field and the Higgs field <£> in terms of six real functions 
a (r,t) , ai(r,i) , a(r,t) , j(r,t) , fx(r,t) and v(r,t) 

1 
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where x is the unit vector in the radial direction. The conditions at the origin are a , a, a\ ■ 
a/r, 7/V, v — > as r — >• 0. The reduced Lagrangian is 
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where the indices run over and 1, ds 2 = dt 2 — dr 2 and we define f^ u = d^a v — d u a^ , 
x = a + 1(7 - 1) , = fi + iv , -D^x = - s)x and = (<9 M - §a M )0. 

This reduced theory has a residual U (1) gauge invariance consisting of gauge transforma- 
tions of the form exp(iQ(r,t)ax/2) with fi(0, t) = 0. Under this, the complex scalar fields 
X and (f> have charges 1 and 1/2 respectively, a M is the gauge field, f^ u is the field strength 
and is the covariant derivative. 

The equations of motion for the reduced theory are 
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The first of these equations is the conserved U(l) current 
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A simplification of the spherical ansatz, hereafter referred to as the reduced spherical ansatz, 
consists of setting the fields a, v and a M to zero, corresponding to parity-even configurations 
($(-x) = $(x), A)(-x) = A)(x), M- x ) = ~A(x))- Thi s leads to the fields 

x = Kv-i) (ii) 

where £ and 77 are real functions of r and t. The equations of motion become 

e_3/n 2 + i_^y' e 

r r \ r ) 2 
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III. SMALL AMPLITUDE EXPANSION 

We look for oscillons in the weakly nonlinear regime, using the weak-field long-wavelength 
expansion. The expansion parameter is e <C 1 with p = er and r = e 2 t as the radial distance 
and the slow time variable respectively. Keeping the scaling general for now, the fields will 
be rescaled as £ = 5^u and i] = 5 v v. The equations of motion for the two fields become after 
rescaling 

d 2 u d 2 u d 2 u d 2 u & 2 

-e 2 5jm 2 H u 3 /2p 2 - e 2 5 2 uv 2 /2p 2 (14) 
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For each region of values of the mass ratio between the two fields (A = <Hir) a different set of 
scaling parametes for the amplitude of the fields can be chosen to lead to oscillon solutions. 
We examine each region separately, following and extending the analysis found in [18] 
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The chosen field scalings in this case are 5$ — 8 V — e. The physical fields are 
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where the complex amplitudes A(p, r) and B(p, r) are defined through the system of coupled 
equations of the Nonlinear Schroedinger type 
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We will call the above set of partial differential equations, envelope equations, since they 
define the evolution of oscillons on large scales and long times. In order to find single 
frequency oscillatory solutions, we use separation of variables with A(p, t) = a{p)e %mHT l 2 
and B(p,r) = b(p)e' tmHT ^ 4: , where a(p) and b(p) are solutions of the following set of ODE's 

^ - m 2 H a + ^b 2 = (20) 
dp 2 Ap 

(21) 
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We will call these profile equations, since they give the time independent (stationary) solution 
to the envelope equations, namely the spatial profile of the oscillon. Finally the physical 
fields are 
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Numerical solutions to the above equations are given in [18] and the oscillons resulting from 
these are found to be extremely long-lived. 

There are more than one ways to treat the case of m H ^ 2m w . One is to take the mass 
ratio to be very close to 2 : 1, that is consider 
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In that case one finds the same form of an oscillon, but the functions a(p) and b(p) are 
given by solutions to a slightly different set of equations. In particular the envelope and 
profile equations for the W field remain unchanged, while the corresponding equations for 
the Higgs beocme 

dA d 2 A 1 P A P A 

2 ' m "^ = -J?-T P B +Am " A (25) 

^ - (m\ + Am%)a + 1& 2 = (26) 

Before concluding this section, we should note the different asymptotic behavior for a positive 
or negative value of Am 2 H . In the former case the two fields decay as a ~ e~ 2mwP /p and 
ft e ~ m wp^ w hil e in the latter case a ~ e ~\/ m H +Am HP and b ~ e ~ mwP , where the in-existence 
of an exponentially decaying solution for Am| < — m 2 H is evident. We will come back to 
this issue later. 

B. mil — 2 * mw = 0(1) 

At this point we have to choose the two scaling parameters 5^ and 5 rj . In [18], both are 

taken to be 0(1), and the resulting fields are to leading order u = Re[A(p,r)e~' tmHt ] and 

v = Re[B(p,T)]e~ imwt } where the complex profiles are defined through the PDE's 

dA d 2 A 3m 2 H \A\ 2 A (6m 2 w - m 2 H )\B\ 2 A 
ZimH dr~ dp 2 2p 2 + 2(m 2 H -Am 2 w )p 2 [ ] 

2im OB _d 2 B _2B m 2 w (Qm 2 w -m 2 H )\A\ 2 B (3m 4 H - 15m 2 H m^)\B\ 2 B 
W dr dp 2 p 2 (m 2 H — Amly) p 2 Am 2 H (m 2 H — Am^p- 

m w( m H - 2m H m w + 2m 2 w )A 2 B 
2mu{jnH — 2mw)p 

where S mHtmw = 1 for m# = mw and zero otherwise. 

The energy of this oscillon is 

E = J drH{r) = ^ J dpH(p) E ~ * (29) 

However this does not seem to be the only choice of scaling parameters. If one chooses instead 
5^ = 0(1) and 5 V = e, one ends up with u = Re[A(p,r)e- imHt ] and v = Re[B(p,r)]e- imwt ] 
where the complex profiles are defined through a different set of PDE's (envelope equations) 

OA d 2 A 14m^ 

- 2lm »8T = W + ^f lAlA (30) 
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dB d 2 B 2B 2mLm 2 



l W" b H I A |2 ; 

dr dp 2 p 2 (16m^m^ — m^) 



-2im w — = — - — + ,_ , w n A M \ B 



+^\A\>B - ±Zrn%\A\*A (31) 

The interesting fact about this scaling choice, is that it leads to a family of solutions, where 
the two fields change in a different way as one changes e. Furthermore, the resulting profile 
equations, when one uses separation of variables in the above set of envelope equations, can 
be solved sequentially as two de-coupled ODE's. Also, choosing different scaling parameters 
for the two fields leads to losing the "Lagrangian structure" of the resulting PDE's, that 
means that they cannot be derived as equations of motion from some effective Lagrangian. 
Finally, the actual physical fields are very different (since one is suppressed by a factor of e). 
Because one of the fields is rescaled by an 0(1) factor, it will dominate the energy, which 
will again be E ~ \ 

We'll call this scaling asymmetric, so as to differentiate between the previous scaling, 
where two fields were scaled by the same number (hence symmetrically). 

The study of the oscillon solutions found in this region will be presented in the Appendix. 



C. A = l + 0(e) 

Finally in this case, using the same general scaling for the field amplitudes, the equations 
of motion become 

^ U l r,2^ U ,4^ U , 2 2^ U , x 2 2o /<-> i ^7 2 /<-> 2r2 2 3 /<-> 2 

di 2 + di&r + ' d7 2 + miiU = t W t5 ^ m H u 3 / 2 P + e J^ V I 9 ~ ' « m " M / 9 

-e 2 5 2 uv 2 /2p 2 + eAm 2 H u + e 2 5^3/2p)Am 2 H u 2 - e 2 5|(l/2p 2 )Am^M 3 (32) 

d 2 v 2 d 2 v A d 2 v 2 2 fd 2 2\ r 2 . 

W + 2e m-r +e dr 2 - + mwV = e {dp-2- 7 2)- + ^m w vu/p 

+e 2 5 v 3v 2 /p 2 - e 2 5 2 v 3 /p 2 - e 2 5^m 2 w vu 2 / p 2 (33) 

Looking at the first equation, we focus our attention on the term eAm 2 H u. This means that, 
whatever the scaling, the 0(e) equation for u will be of the form 

eu htt + m 2 H eui = 0(e) (34) 
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The 0(e) terms are different, based on our choice of scaling parameters 5^ and 5 V . All 
possible such terms are the following: 

5 2 

e5^m 2 H u 2 3/2p + e-V/2p + eAm 2 H u (35) 

To cancel the resonance at this order, the terms on the right-hand side that oscillate with 
the natural frequency of the ODE (ma) have to sum to zero (important note: m# = 2mw 
to zeroth order). Looking at the three terms, the first does not have a part oscillating with 
frequency mn, while the other two do. So these two terms must occur at the same order 
and cancel each other, which leads to the choice 5 2 = 5%. 

The 0(e) part of the expansion now reads (keeping only resonant terms) 

Uijt + m H u i = m 2 H ul3/2p + v 2 /2p + Am 2 H u (36) 
+ m 2 HUl = (l/8p)[2|£?| 2 + (B 2 e- 2imwt + cc)} + {Am 2 H / A){Ae' im « l + cc) (37) 

Canceling the resonant terms gives 

(l/2p)B 2 + (Am 2 H )A = 0^A = -^^B 2 (38) 

If we only wish to use integer powers of e, the first acceptable solution is to take both to be 
0(1). 

The equations hence become 

^ U . r> 2 3 U 4 <9 U 2 2^ U , 2 2o In i 2 In 2 2 3 In 2 

+ 26 ^ + C ^ + = ' ^7 + emHU 3/2/? + ^ /2P ~ 6 /2/? 

-e 2 uv 2 /2p 2 + eAm^w + e 2 (3/2p)Am^ 2 - e 2 (l/2p 2 )Am^w 3 (39) 
d 2 v . „ 2 /;-/• p7 2 /' 2 _ 2 f d2 - \ , 2 



+e 2 3w 2 /p 2 - eV/p 2 - e 2 m 2 y VM 2 /p 2 (40) 

The first order equation for t> is 

+ = 2 m ^ «o/p = (m 2 w /2) [(^e"^-^)' + cc) + ABe"*^^^ + cc)] (41) 

The first term is the resonant one, and canceling it, means setting it to zero, which means 



A = B = 



(42) 
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This means, that in order to avoid resonances at first order, we cannot choose 5^ — 1. The 
next possible integer pair of rescaling parameters is 5^ = e 2 and 5 V = e. These do not create 
any terms for v at 0(e). The equations become in this case 

d 2 u d 2 u d 2 u d 2 u 

W + 2e2 dtd^ + e V + = e V + e3m ^ 23/2p + ey2/2p - e6 < M V2p 2 

-eW/2p 2 + eAm 2 H u + e A (3/2p)Am 2 H u 2 - e 6 (l/2p 2 )Am 2 H u 3 (43) 

d 2 v n 2 d 2 v 4 d 2 v 2 2 f d 2 2 \ „ 2 . 

W + 2t m-r +t dr^ + mwV = 6 W'7J V + e WVU,P 

+e 3 3v 2 /p 2 - eV/p 2 - e A m 2 w vu 2 /p 2 (44) 

From the equation for v, there are no nonlinear terms at order e 2 , meaning that it is im- 
possible to cancel the resonance at this order caused by the term "2e 2 t>/p". This is true 
for all higher order scalings, meaning that canceling the first order resonance in u leads to 
a resonance in v that cannot be canceled. So there is no small amplitude expansion with 
t = e 2 t and p = er and fields rescaled as integer powers of e in this region of the mass ratio. 



IV. NUMERICAL SOLUTION OF OSCILLON PROFILES FOR A = 1 + 0(e 2 ) 

We will now present the solutions to the profile equations 20 and 21 for a range of values 
of the mass mismatch parameter Am|,, obeying the boundary conditions a(p) —> and 
b(p) I p — > as p — > and both decaying exponentially at infinity. 

In [18] the solutions to the ODE's with Am 2 H = were found and fall into a discrete set 
characterized by the number of nodes in b(p). We will focus on the "fundamental" solution 
with zero nodes. Fig. 1 shows the field profiles in the case of the 2:1 mass ratio and it is 
identical to the ones given in [18]. 

We now start to change the mass ratio, by adjusting Am 2 H . It is worth saying that the 
actual mass ratio depends on Am 2 H as well as the expansion parameter e. There are two 
qualitatively different regimes of Am 2 H . On the positive side, there is no a priori upper limit 
for this number. On the negative side, the behavior of the profile equations will change at 
Am 2 H = — 4m ^ =^> A A = — 1, where the linear term in the profile equation for a(p) changes 
sign, and the asymptotic behavior at infinity goes from being exponentially suppressed to 
oscillatory. The existence of this transition in the behavior of the profile equations can be 
explained by the nature of the small amplitude expansion. Let's consider the frequency of 
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FIG. 1: Profiles of a(p)/p (green) and b{p)/p (blue) in the 2:1 mass ratio case. The asymptotic 
behavior of the two curves at p — >• and p — > oo is evident. 




FIG. 2: Profiles of a(p)/p (left) and b(p)/p (right) for different values of the mass ratio. The 
profiles become lower in amplitude and larger in width as the mass mismatch gets smaller, until 
they dissolve for Am^ — > — 1. The change is more dramatic in b(p)/p. 

the Higgs field, both inside the oscillon and in its free state. The squared natural frequency 
of the Higgs field is A4m 2 y = 4m^/(l + AA), while the squared frequency of the Higgs part 
of the oscillon is 4m 2 ¥ (l — e 2 ), which for AA < —1 becomes larger than the frequency of the 
free field. An intuitive necessary condition for the existence of oscillons, regardless of their 
structure or dimensionality, is the requirement that the frequency of the oscillon be lower 
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FIG. 3: Profiles of a{p) (left) and b(p) (right) for different values of the mass ratio. The profiles 
become lower in amplitude and larger in width as the mass mismatch gets smaller, until they 
dissolve for Am 2 H — > — 1. The change is more dramatic in b(p). 

than the natural frequency (or mass in particle terms) of the free field. Thinking of the time 
dependent part of the field as a ball rolling inside a potential well this condition is translated 
into the requirement that the potential gets shallower than quadratic at some point. As long 
as the field can probe this point, oscillons can exist, although their stability requires further 
analysis. It thus follows that if AA < — 1 oscillons cannot exist, which is automatically 
captured by our envelope equations, since they do not admit localized solutions in that 
regime. 

We explored the space — 0.9 2 < Am 2 H < 4 2 , where we were able to find well behaved 
solutions, at least for the fundamental zero-mode solution. The behavior of profiles of the 
two fields are shown in Fig. 2 and 3. Qualitatively, the profiles become larger as the mass 
mismatch parameter gets larger. Furthermore b(p) is more affected by the change in the 
mass than a(p), although the ODE for a(p) is the one that explicitly depends on Am 2 H . 
The results seem to indicate that there is no obvious upper cutoff on the positive side of 
Am 2 H and that we can approach arbitrarily close to the negative cutoff Am 2 H = —Am^y. 
The profiles for both fields become broader and their amplitude decreases as we lower the 
mass mismatch parameter Am 2 H . That was expected, since the profiles seize existing at 
AA = —1, so we are seeing a smooth transition from localization to derealization. Since in 
our numerical code we require the fields to decay exponentially at infinity, derealization is 
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FIG. 4: 1-node profiles of a{p)/p and b(p/p) for Am| = — 0.5 2 ,0, 1. The apparent relative sign 
between the profiles of the two fields is of no physical importance, since the Higgs oscillates with 
twice the W frequency, hence the relative sign of the amplitudes changes twice every W period. 



manifested in the form of a vanishing amplitude. 

In order to see whether or not the behavior of the solutions for changing Am 2 H is accidental 
for the fundamental solution or generic, we solve the profile equations for the single-node 
case. Although we didn't explore the same region of Am 2 H , the existence of solutions and 
their qualitative behavior is the same as before, as is shown in Fig. 4. In this case, the 
solutions continue to grow larger for increasing Am^ and b(p) continues to be affected more 
by the value of the mass mismatch parameter than a(p). 

As we move along the A line for a fixed e the initial energy of the oscillon changes, so we 
need to calculate the energy of the oscillon, based on the e-expansion. The energy for the 
reduced set of equations is 



E 



drH 
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4 V r 



(45) 



Only the last term depends explicitly on the mass mismatch A. Furthermore each term has 
a different dependence on e. We integrate each terms separately (based on the numerical 
solution obtained for the 2 : 1 case) and keep the e dependence of each terms explicit. The 
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reason for using the 2 : 1 profile is that it makes comparison of our results with [18] easier. 
Before that, let's see the terms that are dominant for small e 
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Taking each field to oscillate with its natural frequency (mass), which is correct to lowest 
order in e and integrating the Hamiltonian to find the energy 
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/ drH = e 
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where r\ = e r] res {p)cos{uiHt) and ^ = e £, res (p)cos(cuwt) 

Evaluating the two integrals numerically and neglecting a constant factor multiplying the 
whole energy, but concentrating on the relative size of the two terms, we get 



E ~ e(l + 1.15A) 



(48) 



From this it is clear how the energy changes as a function of the mass mismatch. Since we 
define the W mass to be constant, the mass mismatch only affects the mass of the Higgs, so 
the above term shows also the distribution of the oscillon energy among the two fields. It is 
important to remember that the interaction energy is higher order in e, and hence does not 
appear in this formula. By calculating the energy of the oscillon and comparing with this 
formula, the relative error is 0(e 2 ) as expected. 

V. OSCILLON LIFETIME 



It is now important to evaluate the lifetime of the oscillon for the different mass mismatch 
(relative to the 2:1 ratio). In accordance with [18], we will evolve the solution within the 
reduced spherical Ansatz to show stability. If not perturbed, the solution will remain within 
this Ansatz, so this simplification does not affect the results. Obviously the solution could be 
unstable if perturbed outside the reduced spherical Ansatz, but this does not invalidate the 
analysis that only considers this reduced case. In order to check stability in the full spherical 
Ansatz, the small-amplitude expansion will have to be done for the full set of fields. 

The equations that we will be numerically evolving are 



d"d^ + 2\v 2 r 
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d^r)+^(2r)-3r) 2 + r) 3 ) + ^- ~ V = (50) 

The numerical method used is discretizing the above equations using a second order finite 
difference scheme and evolving them explicitly forward in time. We use a finite spatial box 
having a length of several times the oscillon width with a reflecting boundary at infinity. We 
have seen that the total energy inside the box is conserved and no numerical instabilities 
seem to occur. 

Without loss of generality, we choose units such that g = \/2 and v — 1, the same as in 
[18]. The lifetimes will be calculated using the A = 1 profile as initial conditions even in 
the cases of A ^ 1, to be consistent with [18], since explaining the form of the lifetime plot 
observed by Farhi et al is among the scopes of the present paper. 

The lifetime is defined as the time needed for 50% of the initial energy to leave the 
oscillon region (the oscillon region is defined as the region containing 99% of the oscillon 
energy initially). We choose three values of the small amplitude parameter, namely e = 
0.0125, 0.025, 0.05 and the lifetime plot is shown in Fig. 5. It is important to note, that the 
vertical axis is plotted in units of the rescaled "slow" time, r = e 2 t. This graph is identical 
to the one presented in [18] for the e = 0.05 case. 

The fact that the dynamics of the oscillon within the small amplitude expansion happens 
on the slow timescale r similarly for each e is evident from the lifetime plot. This is a general 
feature of the small amplitude analysis used here, as will be explained later, for all values 
of the mass ratio. Furthermore, lifetimes calculated for three values of e fall almost on the 
same curve, except for a small region of values close to A = 1, which will be analyzed later. 
There one needs to rescale the horizontal axis for each value of e to recover this feature. 

Also we see a deviation of the e = 0.05 curve from the other two for large mass ratio. By 
looking closely into the localized field energy as a function of time for this region, we see 
that the change in lifetime occurs because of a slight change in the decay rate of the Higgs. 
This explains why the discrepancy happens for large mass ratio, where the oscillon contains 
more energy in the Higgs field. In this region a single field Higgs oscillon exists (examined 
in more detail in the Appendix), whose amplitude is e - independent and falls between the 
amplitude of the "regular" two field oscillon for e = 0.025 and e = 0.05. This means that 
the nonlinearity there can be strong enough to affect the decay of the e = 0.05 oscillon, but 
not the other two (e = 0.025 and e = 0.0125). We will not study this slight discrepancy any 
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FIG. 5: Lifetime of the oscillon as a function of the mass ratio A for e = 0.0125, e = 0.025 and 
e = 0.05. Lifetime is measure in units of r = e 2 t and the mass ratio is defined as A = Q "^ H We can 
see that all curves behave identically, with a small deviation of the e = 0.05 curve for large mass 
ratio. The region around ran = 2mjy is not very clear in this graph and will be analyzed later. 

further. 

Perhaps the most striking feature of the lifetime plot is the dip, occurring at a different 
point for each e: for e = 0.05 at A ~ 1.025 and for e = 0.025 at A ~ 1.006. Both values of the 
position of the dip are consistent with A ~ 1 + 10e 2 . This value is taken as an assumption 
for now and will be explained later. 

Apart from the lifetime, it is interesting to see how the energy within the initial oscillon 
region changes with time. We will start with the A > 1 region, since it includes the counter- 
intuitive feature that requires an explanation, namely the dip. We can see two different 
regions in the behavior of the energy evolution, which also occur on the two sides of the dip. 
In Fig. 6 we see two qualitatively different regimes, each of them also occurring on either 
side of the dip. For A > 1 + 10e 2 = 1.025 the decay curves are almost the same, but with a 
shifted origin. 

For values of 1 < A < 1 + 10e 2 , the decay curves are qualitatively different from each 
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other, having a greater decay rate as A increases. On the other hand, when A leaves this 
region (which coincides with the location of the dip), the behavior changes. The initial time 
of constant localized energy is the same, and so is (almost) the rate of energy loss. The 
only difference is the fact that oscillons with different value of A and the same profile have 
different energies. So it simply takes more time to radiate more energy with the same rate, 
that is why oscillons at larger A seem to live longer. This is the observational explanation 
of the dip, the actual physical reason for this behavior will be analyzed in the next section. 

Looking at the region of A < 1, the localized energy as a function of time is shown in 
Fig. 7. For A < 1 — 10e 2 = 0.975 the decay curves are identical, except from being shifted 
by the different value of A. The plateau of the localized energy and the energy loss rate are 
identical for all mass mismatches. As we move into the A = 1 — 0(e 2 ) region, the plateau 
remains the same. However the energy loss rate reduces as we approach the stability point 
A = 1. The oscillon for A = 0.998 = 1 — 0.8 * e 2 shows only an initial loss of energy. This 
is due to the fact that we start with the oscillon profile of the A = 1 case, which has more 
energy than the actual oscillon for A = 0.998. What we see is the initial lump shedding its 
excess energy and relaxing to a stable oscillon configuration. 

In Fig. 8 we compare the decay curves of two oscillons outside the A = 1 + 0(e 2 ) region 
in both sides of the stability point. We see that, once shifted, they practically fall on top of 
each other! 

VI. UNDERSTANDING THE DIFFERENT REGIONS OF THE e — A PLANE 

Now that we have a better understanding of the lifetime plot presented initially in [18], we 
can analyze the whole parameter space for different oscillon behavior. We will use numerical 
simulations both to guide our intuition and to verify our analytical arguments. 

A. Large e 

The analysis of the oscillons is based on the fact that the field oscillates back and forth 
in an anharmonic potential well. However, for large enough value of the field, it will cross 
over to the second part of the double-well Higgs potential. Between the two fields, the W 
(in our small amplitude analysis, the x) wm cross over first. The condition for the field to 
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FIG. 6: Localized energy as a function of time for e = 0.05 and a heavy Higgs. From top to bottom: 
A = 1.5,1.3,1.1,1.05 (blue) and A = 1.001,1.002,1.005,1.015 (red). Oscillons with a large mass 
ratio lose energy at the same rate and larger mass ratio oscillons have a larger initial energy. For 
oscillons close to the stability point (run = 2mw) the energy loss ratio goes down as one moves 
towards the 2 : 1 mass ratio. 

be confined within one well is 

e max — FT/ yi (51) 

max [o(p)\ 

For the oscillon at the 2 : 1 mass ratio, the maximum b(p) is 13.63, so the maximum e 
is 0.073, which is larger (but not by much) to the value 0.05 that was initially used in our 
simulations as well as in [18]. We have seen that b(p) grows with Am 2 H , so e max will decrease, 
as shown is Fig. 9. 

B. |A — 1| > 0(e 2 ) region 

We have seen that lifetime plots include a dip at around A = 1 + 10e 2 . Moreover, beyond 
that point the lifetime plots have a characteristic behavior. There is a plateau until t ~ 800 
or r « 2. Most importantly the decay of the localized energy is almost identical in all 
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FIG. 7: Localized energy as a function of time for e = 0.05 and a heavy Higgs. From top to bottom 
A = 0.998, 0.997, 0.99, 0.98, 0.95, 0.9, 0.8, 0.7. Again the decay rate is growing as one moves 
away from the stability point and freezes, once we exit the A = 1 + (e 2 )region. From that point 
onward the decay curves are shifted copies of each other. 

cases, for values of A on both sides of the stability point A = 1. 

Let's think of the two fields that comprise the oscillon, in the spirit of Feynman diagrams. 
Since we have restricted the fields to be in the reduced spherical Ansatz, their interaction 
has been modified, so we can't easily do actual calculations, but we can use the full SU(2) 
Feynman diagrams to do an order of magnitude estimation. 

By examining a Higgs in its rest frame, assuming it has more than twice the mass of the 
W, the energy-momentum conservation gives us 



where k is the wavenumber of the outgoing W's and c = 1. 

Quantum mechanically the above calculation will immediately give us the lifetime of the 
Higgs. However following the reasoning of [6], we realize that in classical field theory, a 
decaying Higgs field will excite the appropriate wavenumber of the W-field. Hence, these 
wavenumbers have to be present in the original configuration. 




(52) 
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FIG. 8: Comparison between decay curves for A = 0.8 (black) and A = 1.3 (red) and ,A = 1.5 (blue) 
shifted to a common origin. We see that once shifted the three curves are almost identical. All 
three belong to the A = 1 + 0(1) region and their similarity indicates a common decay mechanism 
that is independent of the Higgs mass. 

In the small amplitude analysis, the oscillon has a width of 0(1/ e), which means that 
its Fourier transform has a width of 0(e). By assuming there are no wavenumbers beyond 
some 0(e) value (actually they are there, but are exponentially suppressed), we get 

^ = A-l^A = l + 0( e 2 ) (53) 

So, in order to get an effective pumping of the W-field, we need to be in the A = 1 + 0(e 2 ) 
region. Outside this, even a super-heavy Higgs will classically not decay into W's. Even 
more, the more massive the Higgs, the more exponentially suppressed its decay will be! 

This does not apply only to the heavy-Higgs case. Turning the argument around, one 
can construct a Feynman diagramm including a 4-W vertex, in addition to the Higgs-W 
interaction. The Feynman graph will include three ingoing W's and an outgoing Higgs and 
an outgoing W. Hence, classically we will see a pumping of the Higgs field with a parallel 
demise of the W field. If we take the initial particles to be at rest, energy-momentum 
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FIG. 9: Maximum e for the expansion to be valid, as a function of Ara|. Oscillons that start 
above this curve will have a large enough amplitude to probe both minima of the double well Higgs 
potential. 

conservation can be written as 
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(54) 
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For the long wavelength modes that are present in the oscillon we can Taylor expand the 
above solution for k m w 1, which gives 
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2m 



w 



1 - 0(e 2 ) 



(55) 



As before, the field interaction is restricted to the region of masses with A = l + 0(e 2 ). Since 
the two fields do not interact for |A — 1| > 0(e 2 ), each field must decay on its own, due to 
dispersion. For a gaussian wavepacket in a massive field theory we can write down a general 
expression for the evolution of the width a(t) that is 

t 2 



(56) 



mcr(0) 2 

For a general wavepacket, the fact that the dispersion will be inversely proportional to the 
mass and the square of the initial width will persist. In our case this means that 
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FIG. 10: Evolution curves for A = 1.5 (left) and A = 0.8 (right). We plot the total localized energy 
(blue), the localized Higgs energy (red), the total Higgs energy (purple), the localized W energy 
(green) and the total W energy (cyan). The W decays first, roughly twice as fast as the Higgs, 
and the total energy shed is initially due to the W decay. Furthermore the total energy in the two 
fields is separately conserved, pointing towards an extremely suppressed interaction. 

• The W field will decay roughly twice as fast as the Higgs 

• Both will decay on a timescale of -5 which is again the slow timescale r of our small 
amplitude analysis 

Both features can be seen numerically. 

By setting one of the two fields initially to zero and running the simulation, the resulting 
energy-time plots are identical to when we start with both fields, meaning that the two are 
practically non-interacting. Furthermore, we can see the W decaying almost twice as fast 
as the Higgs. The fact that the W decays first and that the mass of the W is fixed, means 
that initially all decay curves will be identical! This explains the observed behavior, that 
when shifted, decay curves for different vales of A fall on top of each other. 

Simply put, the energy of the oscillon changes with A as E ~ 1 + 1.15A where the term 
proportional to A represents the energy in the Higgs field. The W, which has a constant 
energy content for all values of A decays twice as fast as the Higgs, hence dominates the 
decay rate initially. By changing A we put more or less energy in the field that decays slowly, 
meaning that we change the total oscillon energy but not the initial decay rate, so we change 
the lifetime. 
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FIG. 11: Evolution curves for A = 0.985 (left) and A = 1.015 (right) for e = 0.05. We plot the 
total localized energy (blue), the localized Higgs energy (red), the total Higgs energy (purple), the 
localized W energy (green) and the total W energy (cyan). We see a continuous transfer of energy 
between the two fields. A light Higgs (m# < 2mw) is pumped by the W, while a heavy Higgs 
(run > 2m\y) pumps the W. The latter process is also far more efficient, leading to a quicker 
overall decay of the oscillon. 

C. A = 1 + 0(e 2 ) region 

This is the most interesting region, since it is close to the stability point A = 1 and 
furthermore there is a well denned small amplitude expansion. 

Before going into the stability analysis, we can understand the asymmetry around the 
A = 1 point, where we see that in the A > 1 region the lifetime falls much faster than in 
the A < 1 region in Fig. 11. We see two factors that both lead to the very faster demise 
of the oscillon in the heavy Higgs region. For a heavy Higgs (by "heavy Higgs" we mean 
m H > 2mw), the Higgs field decays and pumps the W fields very effectively and then the 
decay of the oscillon happens through dispersion of the W field. In the opposite region 
(mg < 2mw) the Higgs cannot pump the W, but instead the W pumps the Higgs, which 
is a process that happens far less effectively. Then the oscillon decay curve follows the 
dispersion of the Higgs, which is almost twice as heavy as the W, hence disperses with half 
the rate. 

To lowest order in e the fields in this region are £ = eRe[A(p,T)e~ mHt } and r\ = 
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eRe[V(p, r)e~ %mwt \ where the profiles A and B are defined through their envelope equa- 
tions. According to the small amplitude expansion, processes that occur over a distance 
scale of O(-) and over a time scale of O(^) are captured by the envelope equations. Since 
we expect any instability to involve small wavenumbers (that is long wavelengths) and our 
simulations give us a lifetime of O(^), we will use the envelope equation to study the be- 
havior of the oscillon in the A = 1 + 0(e 2 ) region and compare with results from full wave 
simulations for certain values of the expansion parameter e 

Three important observations justify that the envelope equations (that describe the be- 
havior of the oscillon over long time-scales and long wavelengths) is enough to characterize 
the instability around the A = 1 point. 

• The decay for two choices of e, namely e = 0.025 and e = 0.05 occurs over timescales 
of 0(e 2 ), which is exactly the scaling of time in the small-amplitude expansion 

• Thinking about the decay of an oscillon in terms of Feynman graphs, we conclude that 
the important processes in this case involve wavelengths of 0(e), which is exactly the 
scaling of space in the small-amplitude expansion. 

• Zooming in on the A = 1 + 0(e 2 ) region of Fig. 5 and scaling out the e dependence 
from the A axis, the lifetime curves for both e values are similar, as seen in Fig. 12 
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The envelope equations in the region A = 1 + 0(e 2 ) = 1 + e 2 AA = 1 + e 
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We make the substitution A — > Ae tQl ~ and B — > Be lflT / 2 and the envelope equations become 

dA d 2 A 1 

2im H — = -— - —B 2 + Am 2 H A + 2m H {lA (59) 
or op 2 Ap 

2im w d ^- = - d -^- + ^-B- r ^AB + m w nB (60) 
or op 2 p 2 p 

The stationary solutions to these equations (d/dr = 0) are the oscillon amplitudes a(p,Q) 
and b(p,Q). Evem though it seems that we introduced one extra free parameter, in effect 
constructing an extra family of oscillons, we can always rescale the fields and the space time 
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FIG. 12: Lifetime (in r) of the oscillon around the stability point as a function of the mass mismatch 
(Amjj / Airily) for e = 0.025 and e = 0.05. Lifetime is measured in units of r = e 2 t and the mass 
mismatch parameter is Am^, leading to m? H = 4m^ + e 2 Am|^. We can see that once rescaled the 
curves follow each other very closely, since both follow the same set of envelope equations. 

variables by a constant factor, in order to set Q = rnw- In order to apply the stability 
results as they are found in the literature, we will keep Q explicit for now. 

Following the reasoning found in [24] and [25], which is an extension of the Vakhitov- 
Kolokolov criterion for soliton stability [26] and [27] , oscillons are stable with respect to 
perturbations (solutions to the above equation are non-growing) if 

dQs 



> o 



(61) 



where Qs is the particle number invariant 

Q- 



dp{\A\* + \\B\*) 



(62) 



calculated at the stationary solution (A s = a,B s = b). Although this result is obtained 
for a system of PDE's that differs from our own in the fact that ours has space-dependent 
coefficients, we will use the criterion as described above and test its predictions by a full 
wave simulation. We are currently working towards formally extending the validity of the 
stability analysis for a more general class of systems, like the one presented here. 

We calculate the derivative |^ > at the point Q$ — mw and the results are shown in 
Fig. 13. From this plot we can read the value of the critical value of the mass mismatch 
parameter, referred to as the stability point, which is Am 2 H ~ 8. 



26 



40 



5 - 
O 

* -10 - 




-40 



6.5 



7 



7.5 



8 



8.5 



9 



9.5 



10 




FIG. 13: Stability criterion at £l s = raw = 2 ■ ^ e can see a change in sign occuring at 
Antfj ~ 8, signaling the transition from stability to instability. 

Examining the two sides of the 2 : 1 ratio separately and starting from the light Higgs 
side, we have stability for Am 2 H > —Am^ and no oscillon (meaning also instability) for 
Am 2 H < — 4m 2 v . This translates into a stability threshold at A = 1 — e 2 , which for e = 0.05 
is placed at Xthreshoid = 0.9975. This is numerically verified by taking two values very close 
to this point in Fig. 14. The initial conditions are taken to be equal to the oscillon profile for 
A = 1. In the stable region, we see the oscillon shedding some energy initially and settling 
to a stable configuration, while in the unstable region, we see a very quick pumping of the 
Higgs field, which then disperses. 

Doing the same for a heavy Higgs we find the stability threshold at A ~ 1 + 4e 2 , based 
on Fig. 13, which for e = 0.025 is placed at Xthreshoid = 1-0025. This is numerically seen by 
taking two values very close to this point in Fig. 15. The profiles are taken to be the ones 
derived from the envelope equations including the Am 2 H term. For the stable case we see 
that the total energy curve and the localized energy curve for each field and the oscillon are 
indistinguishable from each other. The localized oscillon energy is constant, while the two 
fields periodically exchange energy. For the unstable case, the Higgs transfers all its energy 
into the W, which then disperses. 
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FIG. 14: Evolution curves for A = 1 - 1.2e 2 = 0.997 (left) and A = 1 - 0.8e 2 = 0.998 (right) for 
e = 0.05, the two sides of the A = 1 — e 2 line. Initial profiles chosen to be equal to the profile at 
A = 1. We plot the total localized energy (blue), the localized Higgs energy (red), the total Higgs 
energy (purple), the localized W energy (green) and the total W energy (cyan). We see either a 
pumping of the W field (left) or a small shed of energy followed by a stable oscillon (right). 




FIG. 15: Evolution curves for A = 1 + 3e 2 rj 1.0019 (left) and A = 1 + 4.5e 2 rj 1.0028 (right) for 
e = 0.025, the two sides of the A = 1 + 4e 2 line. Color coding as before. We can see either a stable 
oscillon (left) with a periodic energy exchange among its constituents, or a very efficient pumping 
of the Higgs field followed by a diffusive W decay (right). 
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VII. CONCLUSION 

We provided an extensive investigation of oscillons in the reduced Spherical Ansatz in 
the SU(2) gauged Higgs model and explained their behavior for different values of the Higgs 
mass and the expansion parameter. 

We can divide the e — A plane in two ways, the criteria being either the stability of 
oscillons and oscillon-like configurations or the formal existence of oscillon solutions in the 
small amplitude expansion. 

Initially we divide the regions of the mass mismatch based on the formal existence of 
solutions within the e expansion framework. 

• There is a single small amplitude expansion in the case A = 1 + 0(e 2 ). 

• The region A = 1 + 0(e) does not contain a small- amplitude two-timing expansion. 

• Finally, there are two scalings, leading to two sets of equations in the case ran — 2mw = 
0(1). The symmetric scaling leads to fields with amplitudes large enough to probe 
both minima of the potential. The asymmetric scaling leads to an oscillon consisting 
predominantly by the Higgs. All these formal oscillon solutions are unstable when 
evolved under the full wave equations. 

The region of interest when searching for stable oscillons is the region near the stable 
point, where furthermore there is a well-defined small amplitude expansion, resulting to a 
set of coupled Nonlinear Schroedinger type equations. 

The (e, A) plane can be divided in regions depending on the behavior of the oscillon (or 
oscillon-like) solutions. 

• The large e region defined as } al u , where oscillons are not stable, due to the 

violation of the small amplitude nature of the expansion. 

• The region where the oscillon fields are not interacting, due to the fact that the relevant 
Fourier modes are exponentially suppressed by the form of the oscillon (or any other 
exponentially localized lump of size O(^). All tested oscillon-like configurations in this 
region are found to be unstable, which each component dispersing independently. 
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FIG. 16: The different regions of the A-e plane. The red curve is the e max line, above which 
the small amplitude expansion breaks down. The two outer green curves A = 1 ± 10e 2 define 
approximately the transition from the interacting to the dispersing region. The two inner blue 
curves A = 1 — e 2 and A = 1 + 4e 2 define the boundaries of the oscillon stability region. 

• The region, where there is interaction and possibly a stable phase-locked oscillon con- 
figuration. This in turn is divided into a stable and an unstable part. 

The characteristic dip in the lifetime graph of the oscillons, corresponds to a change in 
the behavior of the oscillon at A ~ 1 + 10e 2 . This is the point where the A = 1 + 0(e 2 ) 
approximation starts to lose its validity. 

An overview of the different regions is found in Fig. 16 

We are currently applying our insight of this system into simpler models of composite 
oscillons, where the resulting equations have a simpler structure, leaving more room for 
exact Feynman graph calculations and a more rigorous stability analysis. 
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Appendix A: Study of oscillons for mu — 2m\v = 0(1) 
1. Symmetric scaling 

In the case when the mass mismatch is of order one, one uses the appropriate scaling 
analysis leading to the rather more complicated equations described in the previous section. 

Starting with the symmetric scaling, we only numerically computed the node-less (fun- 
damental) solution in the region of A varying ftom 1.5 to 8. With the value of mw = 
this leads to 0.3 ;$ mu — 2mw ~ 2.6. The resulting profiles are plotted in Fig. 17 and Fig. 
18. We see that in both ways of plotting b(p), the curves are smoothly changed as A varies. 
On the other hand, in the case of a(p) such a smooth change is not observed, especially 
when one plots a(p)/p. 

The amplitudes for all values of the Higgs mass in this region exceed unity and are not 
rescaled by a small amplitude parameter. This leads to the fields probing both minima of 
the Higgs potential, leading to the quick decay of the oscillon. 

2. Asymmetric scaling 

For the asymmetric scaling, we again calculated only the node-less solution for mu = 
1,2,3,4, shown in Fig. 19. As the mass increases both field amplitudes get smaller. The 
amplitude in the Higgs field is larger, and when taking into account that the Higgs amplitude 
is not e-suppressed, this is almost a single-field oscillon. 

In the mu — 2mw = 0(1) region the asymmetric scaling 5^ = 0(1) and 5 V = e is only 
the first in an infinite ladder of scalings, namely 5^ = 0(1) and 5 V = e n where n > 1. The 
resulting equations from all these scalings are the same, which means the resulting physical 
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FIG. 17: Profiles of a{p) (left) and a{p)/ p (right) for different values of the mass ratio, 1.5 < A < 8. 
Color-coding goes as follows: 1.5 red, 2, 2.5, 3 blue, 4 green, 5 black 6 yellow 7 cyan 8 magenta 




FIG. 18: Profiles of b(p) (left) and b{p)/ p (right) for different values of the mass ratio, 1.5 < A < 8 

oscillon is the same for the Higgs field and the same rescaled by a different power of e for the 
W field. All these scalings point to an oscillon composed predominately by the Higgs field. 
We will therefore examine the limiting case of this ladder of scalings, n — > oo. Simply put, 
we will set the W field to zero and consider the Higgs in its double-well potential. The only 
difference between this oscillon and other single scalar field oscillons in double well potentials 
found in the literature is that in our case it has space-dependent coefficients, resulting from 
the dimensional reduction due to the use of the Reduced Spherical Ansatz. 
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FIG. 19: Profiles of a(p) (left) and b(p) (right) for m H = 1,2,3,4. 



The equation of motion of the Higgs field in the absence of the W field is 



d^dJ + 2Xv 2 r 



r 2 \r ) 2 [r 



(Al) 



This can be derived from the equation of motion for the physical Higgs field 
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by the substitution 0=^§(l~r) 

Following the e-expansion the £ = Re{A(p, r)e~ imHt } where A(p, r) is defined by 
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By using separation of variables A(p, r) = a(p)e the amplitude a(p) is calculated by 
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A similar stability analysis as before, leads to the oscillons being stable for 
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Calculations for various values of tuh and Q show a globally unstable oscillon. 
The physical Higgs field depends on |. By setting 



(A7) 



33 



the equation of motion becomes 

d 2 i d 2 l 2d£ ox 2 \~ 3~ 2 I-, 

This is the equation of motion of a radial oscillon in 3 spatial dimensions with constant 
coefficients. It is well documented that such an oscillon will be unstable. Thus in the re- 
duced spherical Ansatz the Higgs field forms an unstable 3 dimensional oscillon, that can 
be stabilized near the 2 : 1 mass ratio through interaction with the W field. The interest- 
ing question arises, under what conditions an unstable oscillon can be stabilized through 
interaction with a second field. The second field could even have a different dimensionality, 
meaning a radial laplacian with a different coefficient of the first derivative term. There are 
no general results to our knowledge regarding the interaction of two fields with (effectively) 
different dimensionality within an oscillon. Investigations in this direction will be carried 
out in the future. 



= 



(A8) 



Appendix B: Space and Time Scaling 

Even though we used general scaling for the field amplitudes, the scaling for the space 
and time were chosen to be e and e 2 respectively. One of these can be chosen arbitrarily, 
since e does not have any meaning other than being a small number. Let us consider a single 
scalar field with an equation of motion of the form 

d 2 u d 2 u 

— -7T — tt^" + nr^u + {nonlinear terms} = (Bl) 
at 1 or 1 

If one chooses p = er, the e.o.m. becomes 

(P'll d 2 u 

-—r — e 2 —^ + m 2 u + {nonlinear terms} = (B2) 
at 1 op 2 

Keeping the scaling of the slow time general, say r = e K t, we get 

d 2 u K du d 2 u 9 d 2 u 9 .„ , 

W + 2e ^ + a^- e a? + mu + - = (B3) 

If we want the slow time and the space dependence to enter at the same order, we must 
choose k = 2, hence r = e 2 t. This conclusion is general and does not change for interacting 
fields and any type of nonlinearity. 
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With this final note we have exhausted all possible oscillons in the reduced spherical 
Ansatz of the SU(2) gauged Higgs model. 
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